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INTRODUCTION 


The  recently  increased  interest  in  space  for  civilian  and  defense 
purposes  has  resulted  in  new  concepts  for  the  hypersonic  vehicles.  Since  most 
of  the  currently  planned  and  future  hypersonic  vehicles  will  be  operating  in 
the  upper  atmosphere,  the  'low  density  effects'  will  play  a  major  role  in 
establishing  the  lift,  drag,  moments,  and  aerodynamic  heating  of  a  hypersonic 
vehicle.  In  phase  I  of  this  study,  these  low  desnity  effects  were  analyzed 
for  axially  symmetric  flows.  These  effects,  generally  referred  to  as  the 
surface-  and  shock-slip  boundary  conditions,  were  incorporated  into  a  newly 
developed  viscous  shock-layer  perfect-gas  code.  The  results  obtained  from 
this  code  were  quite  accurate  and  compared  well  with  the  experimental  data. 
Phase  II  work  is  intended  to  analyze  the  low-density  flow  effects  for  bodies 
at  angle  of  attack  with  finite-rate  chemistry  in  the  flowfield. 

Tree,  et  al.,^  included  the  low-density  effects  in  analyzing  the 
hypersonic  ionizing  viscous  shock-layer  flow  past  axially  symmetric  bodies. 
This  analysis  was  limited  to  an  equilibrium  catalytic  wall  and  contained  the 
errors  of  reference  2  (as  explained  in  reference  3)  in  the  specification  of 
surface  slip  boundary  conditions.  The  shock-slip  boundary  condition  also  did 
not  account  tor  the  derivatives  of  the  shock  quantities  in  the  shock-oriented 
coordinate  system  among  other  errors.  This  introduces  errors  in  analyzing 
flows  past  slender  bodies.  Swaminathan,  et  al.,^  and  Song,  et  al.,^  have 
considered  the  surface-  and  shock-slip  effects  for  the  three-dimensional 
flows.  Their  shock-slip  formulation  contains  the  error  mentioned  earlier, 


whereas  their  surface-slip  model 

is  taken 

from 

reference 

6, 

which  also 

contains  errors  as  explained  in 

reference 

3. 

Further , 

the 

frozen-flow 

approximation  employed  in  the  shock-slip  boundary  conditions  of  references  4 
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and  5  i  .at  a  very  physically  appealing  approximation  for  a  thick  shock  under 
the  low  density  conditions,  where  the  characteristic  flow  times  are  large 
enough  to  initiate  chemical  activity  through  the  shock.  Unlike  references  5 
and  b,  it  is  also  desirable  to  solve  the  flowfield  equations  with  equal  step 
sizes  in  the  computational  plane  to  keep  the  same  order  of  accuracy  for  the 
computed  results.  Finally,  to  make  the  computer  code  robust  and  to  speed  up 
global  convergence,  the  input  shock  shape  for  the  low-density  nonequilibrium 
viscous  shock-layer  calculations  should  be  obtained  from  a  corresponding 
approximate  flowfield  calculation  in  place  of  using  an  equilibrium  shock 
shape. Thus,  there  is  a  need  to  develop  an  accurate  and  robust  viscous 
shock-layer  code  for  the  three-dimensional  flows  with  nonequilibrium  chemistry 
and  physically  realistic  surface  and  shock-slip  boundary  conditions  for  a  wide 
range  of  flow  conditions. 

Present  work  was  undertaken  to  study  the  low-density,  nonequilibrium 
three-dimensional  flowfields  around  axially-symmetric  bodies  at  angle  of 
attack.  This  study  was  carried  out  to  develop  a  viscous  shock-layer  code  with 
the  recently  obtained  surface-  and  shock-slip  boundary  conditions.  The 
results  are  obtained  with  an  analytic  grid  built  into  the  equations  with  zero 
numerical  dissipation  of  its  own.  Latest  transport  properties  are  employed 
for  the  air  species  with  appropriate  mixing  laws  for  the  ionized  mixture. 
Detailed  comparison  is  made  with  other  numerical  predictions  (including  the 
Direct  Simulation  rlonte -Carlo  results)  and  the  available  flight  data.  There 
are  not  enough  data  available  for  the  three-dimensional  nonequilibrium  flow, 
especially  under  the  low-density  conditions. 
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ANALYSIS 


Since  its  first  application  by  Davis,^  the  viscous-shock-iuyer  (VSL)  flow 
model  has  been  extensively  used  for  numerical  solution  of  hypersonic  flows 
over  bodies  for  a  range  of  body  shapes  and  freestream  conditions.  Because  of 
the  simplicity  of  VSL  technique  with  the  capability  of  modeling  fairly 
detailed  flowfield  physics,  many  physical  models  have  been  incorporated  into 
VSL  codes  for  a  variety  of  flow  problems  and  detailed  flowfield 
calculations.  Applications  of  VSL  codes  include  predictions  of  flowfields 
around  planetary  entry  bodies  with  massive  ablation  and  radiation,^ 
transitional  and  turbulent  flows  over  slender  vehicles, and  flows  with 
nonequilibrium  chemistry  around  complex  reentry  vehicles.^  The  VSL  technique 
is  quite  attractive  because  of  the  computational  advantage  it  offers.  The 
viscous  shock-layer  equations  are  obtained  from  the  steady-state  Navier-Stok.es 
equations  by  keeping  terms  up  to  second-order  in  the  inverse  square  root  of 
the  Reynolds  number  e .  Consequently,  one  set  of  equations  is  solved  for  both 
the  inviscid  and  viscous  regions.  Further,  parabolic  approximations  are  made 
in  the  strearawise  direction.  In  the  case  of  a  three-dimensional  flow,  the 
crossflow  momentum  equation  is  also  parabolized.  Thus,  the  equations  can  be 
solved  in  both  the  streamwise  and  crossflow  directions  by  marching  techniques 
which  are  efficient  in  terms  of  computer  time  and  storage  requirements.  The 
equations  are  limited,  however,  to  attached  flows,  in  both  streamwise  and 
crossflow  directions,  unlike  some  higher  order  approximations. 

A  3D  nonequilibrium  VSL  code  was  originally  developed  in  reference  10. 
Subsequent  sfudy^  indicated  the  need  for  modifications  and  corrections.  The 
modified  version  of  reference  11  serves  as  the  basis  for  this  study. 
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Governing  Equations 


The  viscous  shock-layer  equations  for  a  reacting  multi-component  gas 

1  O  TO 

mixture  are  obtained  from  the  basic  conservation  equations  >  by  the  method 
mentioned  earlier.  These  equations  in  the  orthogonal,  body-oriented 
transformed  coordinates  (fig.  1)  and  normalized  form  are  provided  here.  Most 
of  the  flowfield  variables  are  normalized  with  their  local  shock  values  to 
keep  them  of  the  order  of  one.  The  normal  coordinate  is  normalized  with  the 
local  shock  stand-off  distance  to  obtain  a  constant  number  of  grid  points 
between  the  body  and  shock.  The  second-order  partial  differential  equations, 
lamely,  5  (or  s)-momentum,  5  (or  $  )-momentum,  energy,  and  the  species 
continuity  equations  are  written  in  the  following  form: 

32w  (ao  3  2§/3n2  +  al  ag/3n)  aW  a2  a3 

3n2  A0(3g/3n)2  3n  AQ(3g/3n)2  AQ(3g/3n)2 

+  — ^-23dr  +  — 9^  * 0  (1) 

A0Og/3n )  H  A0(3g/3n)Z^ 

-  2—2 

where  (3g/3n)  and  3  g/3n  are  the  first  and  second  derivatives  of  the 
stretching  function  g(n).  The  quantity  W  represents  u  in  the  £ -momentum 
equation,  w  in  the  £ -momentum  equation,  T  in  the  energy  equation,  and  in 
tne  species  continuity  equation.  The  coefficients  Aq  to  A^  are  written  as 
follows : 
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nlh  h3  3n  3n 


2  -  -2 


ushp  U  3hl  .  psh  ushp  u  -3nsh  3hl 


hl  h3 


Psh  n  3  °sh  3  p  3  g 
nsh  h3  3^  3^~  3^ 


nsh  hl  h3  n  3p  3n  3^  h3 


n sh  hi  h3  3r| 


3  p  Psh 


A.  =  p  ,  u  .  p  u  /h. 
4  ^  sh  sh  1 


A.  =  p  p  w/h 
5  sh  3 


Energy,  W  =  T: 


2  K.h  = 

e  K 
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A1  - 


p.pvCC  p  u  u  u  C  ,puCn  3o 

sh _ p  p ,  sh  _  sh  sh  ptsh  M _ p  sh 


sh 


h,  n 
1  sh 


35 


p.pwC  ,  C  n  3n  ,  „  K  „  ^  K.  ,  3h., 

M  sh  p _ p,sh  p  sh  _  2  sn  3K.  3_g.  _  2  sh  —  _ 3 

h_  n  ,  3?  ~  e  n  .  3n  -  e  2  . 

3  sh  sh  3n  n  h^  3n 


„  K  K  3  h  2  NS  C  .  w  ,  u  Le,  3C.  3g 
2  sh  _ 1  _  e _  o  p,i  M  sh  1  _ i  _ 


l 

_  i  — 


n2,  h.  3n  nsh  i=l  n  Pr  Pr  3n  3n 
sh  1  sh  sh 


_  NS  w 

A2  ~  W2  "  psh  p  hJ;Tr 


A3 


Psh  Ush  u  3  p  _  Ush  U  P  3  Psh  ^  Psh  Ush  U  11  3  nsh  3  p  3  g 


h,  3  E,  h 


sh 


1  3? 

T  w  n  3 n 


n  ,  h. 
sh  1 


35  3n 


3n 


3_p  3j£  _  W  Psh  3  p  +  Psh  M  n  nsh  3p  3_g 


*sh  3n  3n  '3 
Pa  V 


h3  nsh  3^  3n 


u_u  u  3h  2 

In 


3n 


C.  Oil 

e  2 

"sh 

(  o  ^  oft 

u  ,  r - - 

P  sh  3n  - 
3n 

Oil 

hl 

3n 

2  ush  p 

/3w  3JS.  _  W_ 

3h  2 

• 

e  2 

nsh 

3n  3rT  "  h3 

3n 

+  wL 

A/=Pv,uuPuC  ,  C  /h, 
4  sh  sh  P>sh  p  1 


A=ppwC  C  /h 
5  sh  p,sh  p  3 


With 


NS 

I 

i=l 


h,  w.  =  w.  +  T  w0 
ill  2 


(4b) 


(4c) 


(4d) 


(4e) 


(4f  ) 


(4g) 
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The  remaining  first-order  partial  differential  equations  are  written  as 
follows : 
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Global  Continuity: 


sh  ush  P  u  h3) 


n _  3  nsh  3 _ , 

nsh  35  3n^psh 


u  ,  p  u  ^ — f  h.  h.,  p  p  v)^-^ 

sh  M  3'an  n  an1-  1  3  Msh  M  '  — 

sh  an 


+ 


sT(psh 


P 


w 


hl) 


r\ _  3  nsh  3  ,  — 

nsh  3S  3n  ^Psh  P 


0 


(6) 


n -momentum: 


p  shUstf3  U  a  v  +  p  shp  V  3va_g  _  p  shUshP  U  n  ^  nsh  a  v  3g 

hl  8«  "sh  3“  8^-  '  hl  “sh  3«  3“  8n 

—  —  2 - 2 

+  pshp  w  a_v  pshp  Wp  Snsh  av  ag  psh  uSh  p  u 

h3  ^  “  h3  nsh  3ri  aiT~  hl  nsh  an 

_  0  sh  pw  3  h3  +  |j»h  ip  3_g  =  0 
h3  nsh  an  nsh  3n  a^ 


(7) 


Equation  of  State 

.  . . W 

P  =  P  T  (8) 

W 

It  may  be  mentioned  here  that  the  normal  and  crossflow  velocities  v  and  w, 
respectively,  are  not  normalized.  This  has  been  done  to  avoid  division  by 
small  numbers.  Also,  the  temperature  and  species  concentrations  are  not 
normalized.  The  species  rate  of  production  term  (appearing  in  the  energy  and 
species  continuity  equations  is  written  in  the  forms  given  by  Eqs.  (4g)  and 
(5g)  so  that  the  temperature  and  species  concentration  appear  as  one  of  the 
unknowns.  These  equations  help  in  the  convergence  of  the  numerical  method. 
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The  scale  factors,  h^ ,  and  113 ,  appearing  in  Eqs.  (2)  through  (8)  are 


given  by: 


h^  =  1  +  n 
h2  '  1 

h3  -  rB(l  t,"»sh  cos  9{) 


Further,  Eqs.  (1)  through  (7)  contain  the  following  independent  transformation 


n  -  g(- — •)  -  g(n) 
sh 


(10a) 


where  the  stretching  function  g(q)  is  given  by 


g(n“)  -  1  -  [a  +  - - t°|8. 'n.(yli'l'  I|] 

«. n l (8+1 ) / (6  —  1 ) ]  8  +n(2a+1)  “  1 


(10b) 


and  its  first  and  second  derivatives  are 


dg  _  (  l~q  )(2q+l)  | _ 

<h  anKs+D/^-i)]  [& 


-  n(2a+l)  +  1]  [8  +  n(2a+l)  “  1] 


(10c) 


d2g  _  ( l~a  )  (2q+l)2  ^ _ 

dn~2  £n[(B+l)/(r-l)]  It 


-  n(2a+l)  +  ip  [6+h  (2a+l)  -  1]' 


(I0d) 


Equation  (10b)  permits  the  mesh  to  be  refined  near  the  body  with  o*0  or 
refined  equally  near  both  the  body  and  shock  with  a=  V2  ln  t'ie  physical  plane. 
A  value  of  o  =  V2  has  been  used  here  for  the  low-density  flows.  Parameter  8 
controls  the  amount  of  refinement  with  values  near  1  giving  the  largest  amount 
of  stretching.  Equation  (10b)  may  be  inverted  to  obtain  the  physical 
coordinate  n  from  the  transformed  n  coordinate: 
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(  1~rl~ot'j  ^  iHOLj 

n“-  [1  +  111  (lOe) 

(2a+l)  0-1  0-1 

The  rransf ormat ion  of  tq.  (10a)  ensures  a  uniform  mesh  tr.  the  computational 
coordinate  and  has  zero  built-in  numerical  dissipation  due  to  its  analytic 
nature • 


Boundary  Conditions 

The  following  wall  and  shock  boundary  conditions  are  employed: 


Surface  Slip  Conditions  (employed  at  n  =  n  ) 

-  s 

The  equations  relating  the  slip  values  to  wall  values  and  gradients  at 
the  edge  of  the  Knudsen  layer  (n  =  qg)  are:^ 


Tangential  Velocity  Slip: 


u  =  e  ( 


sh 


^  Pshp  sh  /  p  p 


Psp  s 


(  1 _  3u  8g 

n 


K_  U 


— )1 


W  sh  3n  3n  1  +  n  n  < 
s  sh  £ 


Circumferential  Velocity  Slip: 


(ID 


->  NS  /  W  -  , 

“8  'E  (1/ii/ 


sh 


°°  r  I _ is  _ 


^  Psh P  sh  /  p  p 
s  s 

•)} 


<  w  cos  0 

_ L 


Ws  nsh  3n  3n  1  +  ngh  n  cos  0^ 


(12) 


-11- 


Temperature  Slip  (with  frozen  internal  energy  during  reflection  from  surface): 


sh 


y  -1  s  n  Pr  Pr  ,/p  p  / - 

sh  sh  *sh  Msh  /  PgPs 


sh  Lesh  Les 


i-1  i 


3n 


n  ,  Pr./  p  ,  p  ,  Pr 
sh  sh  sh  sh 


V  NS  W  3/2 

-==^l  (ir)  ct)(lx-Ts  +  IJ 
/  -  -  1=1  1 
/  p  p 

s  s 


NS  3C.  ,  NS  W  8C 

x  y  [  —  _  ( i  _  c  )  y  r  -S-y  —9.  i&i 

A,  L3n  V  i'_4,  3 n 


i-l  ~3n  3n 


q— 1  "q  3n  3rT  s 


Pressure  Slip: 


F5-Fu  +  7^:(^)(f-)  • 


V2, 


Y_1  S  nsh  Prsh  Prs  /psh  psh  '  PS 


1  IT  lei  NS  Wc  3/ 2 

x  1.  Li  <L£.  v  r  — ) 


T  3ri 


tl  (if)  ctIJ 


3n  i=l  i 


sh 


Concentration  Slip: 


cl  ‘  /r  +  **/¥  l 


ush  Lesh,  Le  p 


p  w 

w 


Pr 


sh  Pr  s  /p  ,  p  , 
sh  y  sh 


/  T  p  W  3/2  3C  ( 1  —  C  )  NS  W  3C 

*  /r-V^iirJ  [crsrH  +  -:rr-  i.  rs18*! 


p  ,  p  p  w 
sh  K  s  w  s 


‘sh  8n  an 


,  - .  W  9n  .  —  - 

sh  q-1  q  1  3n  s 
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(13) 


(14) 


(15a) 


where 


3  C .  W.  NS  =  3  C 

c?  =  i  +  (- — -  —  +  rp  n  ip  )  /  y  - aii'j 

1  l3n  3-  Wi  sh^iJs  qiLlWq3n  3-Js 


(15b) 


A  =  0  for  i  =  C>2 

A  =  N  for  i  =  ^ 


^  A  = 


1  /T  Ya  /  p  ,  p  ,Pr  uPr  /  P  P  w  3/2 
—  /  —  f  — ~-A— 'i  /  sh  w  sh  sh  s  /  s  ,  _Sa 

.2/  n  2  -  y/s  u_  ^  -  T  _  -  IwJ 


V  s 


11  s  h  Le  Le  , 
s  sh 


(15c) 


(A  =  0 ,N) 


(15d) 


A  =  0  for  M  =  0^ 
( A  *  N  for  M  * 


(15e) 


(for  all  other  species) 


( 1 5  f  ) 


The  recombination  rate  coefficient,  v  ,  appearing  in  Eq.  (13d)  is  obtained 

A 

from  Scott ^  for  A  »  atomic  nitrogen  (N)  and  from  Zoby,  Gupta,  and  Simmonds15 
for  A  =  atomic  oxygen  (0)  for  finite  catalytic  wall  having  the  shuttle-type 
thermal  protection  system  (TPS).  The  shuttle  TPS  may  be  assumed  noncatalytic 
to  the  recombination  of  other  air  species. 

For  a  noncatalytic  surface,  y  «  0  in  Eq.  (15d)  and  Eq.  (15b)  with 

A 

o 

$ ^  3  0  for  all  the  species  gives: 
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For  a  fully  catalytic  surface,  an  appropriate  boundary  condition  is 


(c.)  =  (C.) 

is  i  <* 


Shock  Slip  Conditions  (at  n  »  1) 


Continuity : 


Osh  vsh  ‘  '  3l“  8 


s-momentum: 


3  n 


^  ysJ^T7  [C°S(B'V  +  ^  sin(g-g  )](^^) 

*  *  *  3n  3n  sh 


sh 


-  (ff^sh  since -0  +  sin  g  ugh  *  cos  g  sin  g 


■u 


where 


and 


®  .h  •  t©Sh  cos<B^)  +  (f^)sh  sin(u-e5)|/Cos  *sh 


(if)5h  -  l(3T^.h  C0S(B V  +  (£)sh  sin<e-S£>)/co,  *sh 


3  U 


sh 


n- mo men turn: 


Psh  =  P»  +  sin  6(sin  g  +  vsh) 
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(17) 

(18) 

(19a) 

(19b) 

(19c) 

(20) 


Energy : 


sin  6 


{“ —  [cos(B-e  )  + 


3n  sh  * 


NS  : 

+  lml  Ci,sh  hi,sh  '  it  rush"  C0S  S)  +  ^Sin2e  '  %h>]  +  iiS" 

u  h  NS  1  3n 

X  ^  ill  L6i’sh  hi,^iriCOS^~Q^  +^sin(6-65)] 

3  C  3  C .  NS 

*  ,ln<6",5)l  ‘  l  ci,.  V 


Species  Continuity: 


e  u^hLe 


fr^fr^V +  Ir1  >](lr 


+  sin  g (C 


,  u  ~  )  *  f  w.dn«wJ 

i,sh  i,oo  J  x  i,ave  sh 

sh  ’ 


(22a) 


w  -  w  /  2 

i,ave  i,sh 


(22b) 


6Sh  *  (nco  -fW 


(22c) 


The  shock  thickness,  {  ,  to  a  good  degree  of  approximation  may  be  obtained 


6  -  0.174  M  e' 

Sh  oo 


(22d) 
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Equation  (22d)  has  been  approximated  from  references  16  and  17.  This  equation 

has  been  further  verified  by  integrating  the  equations  for  a  merged  stagnation 

shock  layer  of  nonequilibrium  dissociating  gas  with  the  same  properties  as 

those  used  here.  The  shock  thickness,  6  ,  ,  obtained  in  reference  18  from 

sh 

integrating  through  the  shock  keeping  the  w_^  term  is  very  close  to  that 

obtained  from  Eq.  (22d).  If  the  shock  is  thin  (i.e.,  6  ,  is  of  order  10-2) 

sh 

the  mass  rate  of  production  of  atoms,  w_^ ,  may  be  assumed  negligible  across  the 
shock.  This  assumption  makes  the  right-hand  side  of  Eq.  (22a)  zero. 

Since  velocity  components  tangent  and  normal  to  the  shock  (see  Fig.  1) 
are  not  the  same  as  those  tangent  and  normal  to  the  body,  the  following 
transformations  are  needed  to  obtain  velocity  components  at  the  shock  in  the 
body-oriented  coordinates  from  the  shock-oriented  coordinates: 


ush  =  “sh  cos  *sh  cos(6~V  '~Vsh  cos  *8h  sin(6_V  <23) 


'sh  =  “sh  cos  *sh  sln(S-V  +  v8h  cos  *sh  cos  (6  ~0£  ) 


(24) 


w  =  -  sin  $ 
sh  sh 


/  ~ 2  ~2 

7  ush  +  vsh 


(25) 


where  £  and  9  are  the  shock  and  body  angles,  respectively,  and  is  the 

19 

crossflow  angle  (see  fig.  2a)  defined  as: 


sin  *  , 
T  sh 


w 

u- 

00 


(26a) 


=  -  sin(o '  -  0')N  +  cos(o' 

L 


92)T, 


(26b) 
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In  Eq.  (26b),  (a'  -  0')  is  the  angle  which  the  velocity  vector  after 

*  A  A 

crossing  the  shock  makes  with  the  shock-tangent  vector,  T,  in  the  N-T  plane 

(see  fig.  2b)  and  N  and  T  are  the  components  of  normal  and  tangential 

4>  4> 

vectors  N  and  T,  respectively,  in  the  crossflow  direction  (see  fig.  2a).  For 

a  2D  flow  over  an  axisymmetric  body  at  zero  degree  angle  of  attack,  Eqs.  (23) 

and  (24)  reduce  to  those  given  in  reference  20  with  d>  =0  degree. 

sh 

Alternately,  the  following  relations  can  also  be  used  to  relate  the 
velocity  components  at  the  shock  in  the  body-  and  shock-oriented  coordinate 
systems : 


u  ,  =  cos 
sh 


I  ,  COS(g  ,  -0 
sh  sh  £ 


,  J~  2  ~2 

)  /  U  ,  +  V  , 

sh  sh 


(27) 


sh 


'  cos  *sh  sln(6sh  ■  e£  “sh  +  ' 


2 

sh 


(28) 


sh 


■  sin*si/“L  + 


sh 


(29) 


where  B  is  the  flow  angle  from  the  body  axis  in  a  meridional  (or  41  * 
constant)  plane^^  (see  fig.  2c): 


-sin(a'  -  0')N  +  cos(a'  -  9')T 

tan  R  =  — —  - - — 1 - _ £ 

sh  -sin(o  -  0  „  )N  +  cos(a  -  e')T 

»  X  Z  X 


(30) 


With  angles  a'  and©'  defined  in  figure  2b  and  Nx,  Tx,  and  Nr,Tr  are  the 
components  of  vectors  N  and  T  in  the  body-axis  and  radial  directions, 
respectively . 

Angle  (8sh“0^)  appearing  in  Eqs.  (27)  and  (28)  may  be  obtained  from  (see 
fig.  2c): 


ta»(Bsh  -  0?) 


sh 

j* 

sh 


sh 

*sh 


(31) 
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Chemical  Reaction  Model 


The  net  mass  rate  of  production  of  a  chemical  species  per  unit 

*  9  1  9  9 

volume,  w*,  is  obtained  from  the  usual  kinetic  relations.  »  For  the  seven 
species  (C^ , ^ ,0 ,N ,N0 , N0+,  and  e  )  air  calculations  in  the  present  code,  the 
reaction  model  is  similar  to  that  used  in  references  22  and  23.  The  species 
production  terms  appearing  in  the  energy  and  species  continuity  equations  are 
written  in  terms  of  temperature  and  concentration  following  the  approach  of 
references  21  and  22.  The  reaction  rate  constants  and  the  relative 
efficiencies  of  the  catalytic  third  bodies  are  obtained  from  reference  22. 
These  rate  equations  and  constants  are  similar  to  those  used  in  previous 

studies  These  are  also  identical  to  those  used  in  2D  VSL  studies 

15,24_ 

Thermodynamic  and  Transport  Properties 

Thermodynamic  properties  for  the  multicomponent  air  mixture  are  obtained 

by  assuming  a  thermally  perfect  gas.  The  equation  of  state  for  the  gas 

mixture  is  given  earlier  as  Eq.  (8).  Thermodynamic  properties  for  the 

2S  27 

individual  species  are  obtained  by  interpolation  in  the  data  of  Browne.  ’ 
More  details  of  the  formulation  employed  in  the  present  calculations  can  be 
found  in  reference  23,  which  follows  the  thermodynamic  models  of  references  25 
to  27  . 

Transport  properties  for  the  gas  mixture  are  obtained  using  the  methods 
of  reference  28  for  viscosity  and  that  of  reference  29  for  thermal 
conductivity.  Individual  species  viscosity  and  frozen  thermal  conductivity 
are  based  on  the  more  recent  work  of  reference  30.  These  data  are  believed  to 
be  more  accurate  at  higher  temperatures,  encountered  in  nonequilibrium 
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calculations.  It  may  also  be  mentioned  that  the  present  mixture  thermal 
conductivity  relation^  accounts  for  the  contribution  from  the  internal  modes 
similar  to  reference  11,  whereas,  the  relation-^  used  in  reference  31  does 
not . 

In  this  study,  a  variable  Prandtl  number  is  used,  whereas,  the  Lewis 
number  is  set  equal  to  1.4.  For  the  ionized  species,  ambipolar  type  of 

O') 

diffusion  is  employed.  For  such  a  diffusion  the  binary  diffusion 

coefficients  must  be  doubled  in  the  calculations.  Accordingly,  a  Lewis  number 
of  2.4  has  been  used  for  the  ionized  species.  Previous  studies^ ignored 
this  effect. 

Method  of  Solution 

The  numerical  solution  of  the  viscous-shock-layer  equations  was  presented 

2  21 

by  Davis^*^  for  the  axi symmetric  and  two-dimensional  (2D)  flows.  An  implicit 
finite-difference  scheme  was  used  which  was  efficient  in  terms  of  both 
computer  time  and  storage.  In  developing  the  three-dimensional  (3D)  viscous- 
shock-layer  method  for  a  perfect  gas,  Murray  and  Lewis^  adopted  a  similar 
scheme  for  the  numerical  solution,  as  did  references  4,5,11,  and  31  for  the 
nonequilibrium  three-dimensional  formulation.  The  same  finite-difference 
method  is  used  here  to  solve  the  5 -momentum,  £ -momentum,  energy,  and  species 
continuity  equations,  namely,  equations  (1)  through  (5).  By  using  two-point 
backward  differences  in  the  5 -direction,  the  zig-zag  finite-difference 
expressions  of  Krause^  in  the  5 -direction,  and  central  differences  in  the 
n -direction  in  the  second-order  parabolic  equaitons  (1)  through  (5),  a 
tridiagonal  finite-difference  form  is  obtained  which  can  be  solved  by  the 
Thomas  algorithm  for  each  of  the  dependent  variable  represented  by  W.  The 
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solution  to  this  second-order  set  of  equations  is  obtained  in  an  uncoupled 
cascading  scheme  similar  to  that  used  by  Davis. ^  The  first-order  differential 
equations  (6)  and  (7)  of  continuity  and  normal  momentum,  respectively,  are 
solved  by  a  similar  method,  but  they  are  coupled  together. ^  The  coupled 
solution  of  these  two  equations  promotes  converence  and  ensures  stability  in 
the  marching  solution,  especially,  for  the  slender  bodies. 

There  are  a  few  points  in  the  flowfield  where  the  governing  equations 
become  singular  and  the  limiting  form  must  be  evaluated.  Along  the  stagnation 
streamline,  the  5 -derivatives  of  the  pressure  and  shock,  stand-off  distance  are 
zero  and  the  tangential  component  of  velocity  behind  the  shock,  ugh,  is  also 
zero.  These  result  in  a  singularity  in  the  coefficients  of  ^-momentum 

equation.  The  continuity  equation  also  becomes  singular  with  5  approaching 
zero.  There  singularities  are  removed  by  using  the  series  expansions  of 

references  36  and  37.  Further,  the  crossflow  velocity,  w,  is  zero  along  the 
windward  and  leeward  streamlines  and  the  value  of  3w/g 5  is  needed  in  the 
continuity  equation.  This  is  obtained  by  taking  the  <; -derivative  of  the 

crossflow  equation  and  solving  for  the  new  dependent  variable  3w/9f;.  In  the 
windward  and  leeward  planes,  the  ^-derivatives  of  all  dependent  variables 
except  that  of  w  are  zero,  and  the  second  derivative  of  w  is  zero.  The 
equation  for  3w/3£  is  solved  by  writing  it  in  the  standard  parabolic  form, 

namely,  Eq.  (1). 

The  present  viscous  shock-layer  solution  begins  on  the  spherical  nosecap 
at  the  aerodynamic  stagnation  point  where  an  axisymmetric  solution  in  a  wind- 
fixed  coordinate  system  is  obtained.  The  axisymmetric  (one  plane)  solution  on 
the  sphere  is  marched  downstream  to  encompass  the  aerodynamic  and  geometric 
stagnation  points.  The  wind-fixed  solution  is  then  rotated  around  the  sphere 


-20- 


and  Interpolated  to  obtain  the  shock-layer  profiles  in  the  body-fixed  (three- 
dimensional)  coordinate  system.  These  are  used  as  the  initial  profiles  for 
the  3D  solution.  The  3D  solution  begins  in  the  windward  plane  at  a  specified 
streamwise  (or  5 )  location  and  marches  around  the  body  obtaining  a  converged 
solution  at  each  circumferential  (or  ^)-direction  step.  The  governing 
equations  are  solved  at  each  point  in  the  flowfield  in  the  following  order: 
species  continuity,  crossflow  (or  5 )-momentum,  energy,  streamwise-momentum, 
integration  of  global  continuity  for  shock  standoff  distance,  ng^,  and  the 
coupled  normal  (or  q )-momentum  and  continuity  equations.  After  integration  of 
the  continuity  equation,  the  new  shock  standoff  distance  is  used  to  update  the 
shock-layer  profiles  for  the  next  local  iteration  while  the  shock  slope, 
30^/35,  remains  fixed  at  the  input  value.  The  equations  are  iterated  in 
this  manner  until  each  of  the  flowfield  variables  is  converged  through  the 
entire  shock-layer  at  the  given  5-  and  ^-stations. 

The  solution  procedure  at  a  given  5  -station  on  the  axisymmetric  part  (or 
nose  region)  is  shown  in  figure  3.  Figure  4  depicts  the  solution  procedure 
used  on  the  body  (beyond  the  tangency  point)  at  a  given  5 -station  for  various 
values  of  ^  (<£  )  or  cross-flow  planes. 

Global  convergence  of  the  VSL  solutions  is  satisfied  by  repeating  the 
calculation  over  the  entire  body  using  the  computed  shock  standoff  distance 
obtained  during  the  previous  pass  to  provide  the  new  shock  shape  (i.e.,  shock 
standoff  distance  and  its  slope).  Global  iterations  are  repeated  until  the 
input  and  output  shock  shapes  are  converged. 
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Calculation  of  Input  Shock  Shape 


The  global  Iteration  process  outlined  in  the  previous  paragraph  requires 
an  initial  value  of  the  shock  slope,  3  ns C  »  and  3v/9C,  which  introduce  an 
elliptic  effect  into  the  equations.  In  the  present  scheme,  3nsh/H>  Is 
obtained  either  from  an  inviscid  equilibrium  air  solution^®  (which  is  scaled 
by  a  constant)  or  from  a  viscous  nonequilibrium  analysis.  The  values  of  9v/3£ 
are  calculated  from  a  backward  difference  approximation.  The  scaling  factor 
used  with  the  inviscid  equilibrium  solution  for  the  shock  standoff  distance  is 
in  the  range  of  1.2  to  1.25.  The  initial  input  shock  shape  obtained  from  this 
scaling  accounts  for  the  nonequilibrium  effects  on  the  flowfield  and  limits 
the  shock-shape  correction  and,  therefore,  decreases  the  required  global 
iterations. 

Even  though  the  results  obtained  from  using  an  inviscid  equilibrium  shock 
shape  with  a  scaling  factor  are  satisfactory  for  nonequilibrium  viscous  flows 
at  high  Reynolds  numbers,  convergence  problems  are  encountered  at  freestream 
Reynolds  numbers  (based  on  the  nose  radius)  less  than  about  2000.  These 
values  of  the  Reynolds  number  are  encountered  under  the  low-density  flight 
conditions.  Under  these  conditions,  the  entire  shock-layer  becomes  viscous. 
An  alternate  method  for  obtaining  a  3D  viscous,  nonequilibrium  input  shock 
shape  is,  therefore,  suggested  here.  Starting  point  is  the  generation  of  2D 
viscous  nonequilibrium  shock  shapes  for  an  axisymmetric  body  (see  fig.  5)  with 
effective  body  angles  of  (0^-hx)  (for  simulating  the  flowfield  on  the  windward 
side)  and  (0  -tx )  (for  simulating  the  flowfield  on  the  leeward  side).  The  2D 
VSL  solutions  can  be  obtained  from  the  3D  VSL  code  by  running  it  in  the  2D 
mode.  For  the  body  at  an  angle  of  attack,  a>  these  solutions  are  used  from 
the  aerodynamic  stagnation  point  up  to  the  tangency  point  (where  the  nose  of 
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the  body  is  axisymmetric)  in  the  wind-fixed  coordinates  (see  fig.  6).  Beyond 
the  tangency  point,  the  shock  shape  (for  the  body  solutions)  is  specified  in 
the  body  fixed  coordinates.  The  relationship  between  the  wind-fixed  and  body- 
fixed  coordinates  (shown  in  figs.  1  and  5)  along  the  windward  and  leeward 
streamlines,  respectively,  is 


( S*  )  3  C* 

wind's ($ )“0  body 


(32a) 


^Swind^s  (<j>  )“ir  Sbody  +  a 


(32b) 


Having  obtained  the  windward-side  and  leeward-side  values  of  the  shock 
standoff  distance  and  shock  slope  from  the  2D  solutions  of  bodies  with 
effective  body  angles,  their  distributions  around  the  body  (0  <  £  (<j, )  <  xr )  for 
body  solution  (l.e.,  beyond  the  tangency  point)  may  be  obtained  from: 
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The  same  relations  are  valid  for  s*  .  >  (s*  ,  )„„. 

body  '  body  TP 
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Surface  Quantities 


The  following  surface  quantities  may  be  computed  from  the  flowfield 
profiles  obtained  from  solutions  of  the  governing  equations  described  earlier. 


(i)  Stanton  Number: 


St 
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(ii)  Heat  Transfer  Coefficient: 
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(iii)  Wall  Heat-Transfer  Rate: 
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(iv)  Skin-friction  Coefficient  in  the  Streamwise  Direction: 
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(v)  Skin-friction  Coefficient  in  the  Transverse  Direction: 
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(vi)  Pressure  Coefficient: 
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RESULTS  AND  DISCUSSION 

Detailed  flowfield  results  have  been  obtained  here  for  sphere-cone  shaped 
bodies  at  various  angles  of  attack  with  a  3D  nonequilibrium  viscous-shock- 
layer  code  (VSL3DNQ).  Obtained  results  show  the  effect  of  low  density  on  the 
flowfield  around  hypervelocity  vehicles.  Present  predictions  have  been 
compared  with  the  available  experimental  data  and  other  numerical  predictions. 
Including  those  from  the  Direct  Simulation  Monte  Carlo  Method  (DSMC). 

Perfect-gas  Results  With  and  Without  Slip 

Perfect-gas  calculations  were  carried  out  first  to  check  the 
computational  grid  implemented  here.  Figure  6  shows  comparison  of  the 
laminar  heating  distributions  with  the  experimental  data  of  deary  (ref.  39) 
for  a  15-deg  sphere  cone  at  20-deg  angle  of  attack.  The  present  prediction 
compares  quite  well  both  for  the  windward  (?($)“0°)  and  side  (t($)~90°) 
meridior.  ils.  Figure  7  contains  the  surface  heat  transfer  distributions  for 
these  two  meridional  planes  with  and  without  slip.  Similar  to  figure  6  these 
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results  are  also  for  a  15-deg  sphere  cone  at  20-deg  angle  of  attack.  Only 
difference  here  (in  fig.  7)  is  the  value  of  freestream  Reynolds  number,  which 
is  much  lower  due  to  the  low  value  of  the  density  employed  here.  The  Reynolds 
number  parameter,  e ,  which  is  a  measure  of  low-density  effects,  has  a  value  of 
0.138  here.  We  should  mention  here  that  the  axisymmetric  solutions  over  the 
spherical  nose  are  given  in  the  wind-fixed  coordinates  (see  fig-  1),  whereas 
the  3D  (or  body)  solutions  obtained  beyond  the  tangency  point  are  given  in 
body-fixed  coordinates.  Figures  6  and  7  illustrate  that  fairly  accurate  and 
consistent  results  are  obtained  for  the  windward  and  side  meridional  planes 
with  and  without  slip. 

Results  for  Nonequilibrium  Flow  Without  Slip 

These  results  have  been  obtained  to  check  the  accuracy  of  solutions  with 
the  new  grid  and  transport  properties  as  well  as  to  study  the  effect  of  input 
shock  shape.  As  mentioned  earlier,  the  input  shock  shape  plays  a  major  role 
in  obtaining  good  quality  solutions  and  their  global  convergence. 

Rffect  of  input  shock  shape.  Figure  8  shows  the  input  shock  shapes  for 
the  leeward  and  windward  sides  obtained  from  the  inviscid  equilibrium  and 
viscous  nonequilibrium  flowfield  calculations.  As  mentioned  earlier,  the 
inviscid  equilibrium  shock  shape  is  obtained  from  the  code  of  reference  38. 
The  viscous  nonequilibrium  shock  shape  is  obtained  from  2D  calculations  using 
the  effective  body  angles  for  the  windward  and  leeward  sides  and  employing  a 
cosine  distribution  to  obtain  shock  shapes  for  other  meridional  planes.  As 
expected,  the  shock  standoff  distance  obtained  from  the  inviscid  equilibrium 
calculations  is  smaller  than  the  one  obtained  from  the  viscous  nonequilibrium 
computations.  The  difference  between  the  two  values  increases  with  the 
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decrease  in  density.  Figure  9  gives  the  same  shock  standoff  distance  by 

employing  the  two  different  shock  shapes  of  figure  8  for  this  relatively  high 

density  case.  It  may  be  mentioned  here  that  for  moderately  high-density 

nonequilibrium  flowfield  calculations,  an  input  shock  shape  obtained  from 

inviscid  equilibrium  flow  calculations  serves  as  a  good  initial  guess  with  a 

scaling  factor  of  about  1.2.  However,  for  the  fairly  low  density  flows  (R  < 

e  ,oo 

2000  or  so),  solutions  may  never  converge  with  an  inviscid  equilibrium  input 
shock  shape. 

Figures  10  and  11  show  the  convergence  history  of  the  shock  slope  using 
the  viscous  nonequilibrium  and  inviscid  equilibrium  input  shock  shapes, 
respectively.  Clearly,  the  shock  slope  converges  faster  with  a  viscous 
nonequilibrium  input  shock  shape.  Surface  heat  transfer  coefficient  employing 
the  two  input  shock  shapes  converges  to  almost  same  value  as  shown  in  figure 
12,  with  inviscid  equilibrium  input  shock  shape  requiring  an  additional  global 
pass  for  this  case. 

Comparison  with  2D  results.  Figures  13  through  18  show  comparisons  of 
the  results  obtained  for  a  20-deg  sphere  cone  at  zero-deg  angle  of  attack  by 
employing  the  present  3D  nonequilibrium  viscous-shock-layer  (VSL3DNy)  code  and 
an  earlier  developed  (ref.  40)  2D  nonequilibrium  viscous-shock-layer  (VSL2DNQ) 
code.  Both  of  these  codes  employ  similar  computational  grids,  reaction-rate 
kinetics,  transport  and  thermodynamic  properties.  Therefore,  the  two  codes 
predict,  as  expected,  almost  the  same  velocity,  temperature,  and  species 
concentration  profiles  as  shown  In  figures  13  through  18. 

3D  calculations.  Figures  19  through  25  contain  the  3D  results  obtained 
for  a  20-deg  sphere  cone  at  5-deg  angle  of  attack.  Figures  19  and  20  show  the 
skin-friction  coefficient  and  surface  heat-transfer  coefficient  distributions, 
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respectively.  Results  for  the  axisymmetric  spherical  nose  are  shown  in  wind- 
fixed  coordinates  (see  fig.l)  whereas  the  3D  results  for  different  cross  flow 
planes  beyond  the  tangency  point  are  in  body-fixed  coordinates  in  these  two 
figures.  Figures  21  through  25  contain  flowfield  profiles  of  various  species 
at  different  crossflow  (or  meridional)  planes  for  a  given  streamwise  body 
location  of  £ (s)  =  2.02.  For  the  calculations  shown  in  these  figures,  the  wall 
is  assumed  to  be  noncatalytic  and  the  species  are  assumed  frozen  at  the 
freestream  value  behind  the  shock  i.e.,  at  n/ng^=l.  These  figures  also  show 
that  maximum  dissociation  and  ionization  (as  depicted  by  the  electron  number 
density  provided  in  fig.  25)  occurs  along  the  windward  meridional  plane 

(c  (<(>  )  =  °°)  • 

Effect  of  angle  of  attack.  The  shock  slope  and  the  surface  heat  transfer 
coefficient  for  a  9-deg  shpere  cone  at  three  different  angles  of  attack  (5°, 
8°,  and  11°)  are  given  in  figures  26  and  27,  respectively.  The  results  are 
shown  only  for  the  leeward  and  windward  meridional  planes.  As  before,  the 
axisymmetric  results  for  the  spherical  nose  are  given  in  wind-fixed 
coordinates,  whereas  the  3D  results  beyond  the  tangency  point  are  given  in 
body-fixed  coordinates.  The  maximum  heating  at  the  aerodynamic  stagnation 
point  is  obtained  for  the  5-deg  angle-of-attack  case.  Also,  the  heating  is 
higher  both  for  the  windward  and  leeward  sides  of  the  body  at  5-deg  angle  of 
attack  than  at  higher  angles  of  attack. 

Results  for  Nonequilibrium  Flow  with  Slip 

Comparison  with  experimental  data.  Before  carrying  out  the  extensive 
nonequilibrium  calculations  for  the  low  density  flows,  RAM  C-Il  and  RAM  C-III 
(refs.  41  and  42)  flight  data  were  used  to  verify  the  accuracy  of  VSL3DNQ 
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predictions  with  slip.  These  were  two  of  the  three  "Radio  Attenuation 
Measurement"  test  vechicles  for  the  investigation  of  flowfield  plasma  under 
reentry  conditions.  The  vechicle  was  a  9-deg  spherically-blunted  cone  and  was 
instrumented  to  measure  electron  number  density  across  the  shock  layer. 
Microwave  ref lectometers  were  used  to  measure  the  peak,  values  of  electron 
number  densities  in  the  shock  layer  at  several  locations  along  the  cone.  A 
Langmuir  probe  rake  that  extended  across  the  shock  layer  at  the  base  of  the 
cone  measured  elctron  number  density  profiles.  Details  of  the  vehicle  and 
experiment  are  given  in  references  41  and  42.  RAM  C  flight  data  are  the  only 
data  for  the  flowfield  profiles  under  the  high-energy,  low-density  flight 
conditions.  Figure  28  shows  an  excellent  comparison  between  the  flight-data 
and  the  present  predictions  with  surface  and  shock  slip.  Similarly,  the 
present  predictions  of  peak  electron  number  density  within  the  shock  layer 
compare  very  well  with  the  experimental  data  as  shown  in  figure  29.  These 
comparisons  show  that  the  accuracy  of  the  low-density  predictions  with  the 
present  slip-formulation  is  quite  good. 

Comparisons  with  DSMC  and  other  calculations.  Figure  30  gives  z 
comparison  between  the  present  predictions  and  those  of  reference  43  from  the 
direct  Simulation  Monte  Carlo  (DSMC)  method.  Present  results  are  shown  with 
and  without  the  surface  and  shock  slip  boundary  conditions.  These  predictions 
(with  corrections)  are  based  on  the  surface  boundary  conditions  of  reference  3 
and  shock  slip  conditions  presented  here.  These  predictions  also  employ  the 
computational  grid  of  Eq.  (10).  The  results  without  corrections  are  those 
obtained  from  the  code  of  reference  11  with  slip  formulation  of  references  4 
and  5.  These  results  are  used  in  reference  43  for  comparison  with  the  DSMC 
predictions.  Clearly,  the  present  VSL3DNQ-slip  predictions  (with  corrections) 
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are  in  fair  agreement  with  the  L)SMC  results  and  have  the  similar  behavior  at 
high  altitudes.  The  differences  between  two  results  may  be  due  to  the  surface 
interaction  models,  the  transport  properties,  and  the  computational  grid 
employed  in  the  calculations.  Present  predictions  as  well  as  the  DSMC  results 
approach  the  no-slip  VSL  predictions  under  the  high-density  flow  conditions, 
as  expected.  We  should  mention  here  that  substantial  computer  resources  (both 
computational  time  and  storage)  are  required  to  obtain  the  DSMC  results, 
especially,  at  moderately  high  densities  such  as  those  shown  between  60-  and 
75-km  altitudes  in  figure  30.  Comparison  between  the  present  predictions  and 
DSMC  results  (of  ref.  43)  for  the  surface  pressure  coefficient  is  given  in 
figure  31.  Unfortunately,  only  one  computed  result  is  available  from  the  DSMC 
method  for  the  conditions  of  this  figure.  Still,  the  predictions  compare 
favorably . 

Figures  32  through  34  show  further  comparisons  with  the  code  of  reference 
LI  (containing  the  slip  formulations  of  refs.  4  and  5)  and  the  present 
results.  Present  results  are  shown  'with  corrections'.  Present  slip 
formulation  gives  lower  flowfield  temperatures  and  lower  surface  heat-transper 
races,  particularly,  in  the  stagnation  region. 

Low  density  effects  on  3D  calculations.  Finally,  detailed  flowfield 
results  are  obtained  to  analyze  the  effects  of  low  density  on  3D  calculations. 
Figures  33  through  46  contain  these  results  for  a  9-deg  sphere-cone  at  5-deg 
angle  of  attack  for  a  fully  catalytic  surface.  From  the  results  presented  in 
figures  33  through  38,  the  surface  pressure  coefficient  (fig.  38)  is  impacted 
the  least  by  the  slip  effects.  The  effect  of  slip  on  other  surface  quantities 
is  also  not  much  due  to  the  fully  catalytic  wall  boundary  condition  employed 
here.  Of  all  the  flowfield  quantities,  temperature  (fig.  40)  and  density 


-30- 


(fig.  41)  appear  to  contain  the  maximum  slip  effects  near  the  surface. 
Species  mass  fractions  near  the  fully  catalytic  surface  are  the  same  as  those 
in  the  freestream  as  shown  in  figures  42  through  46. 


CONCLUDING  REMARKS 

Results  from  this  study  can  be  summarized  as  follows: 

1.  A  3D  nonequilibrium  viscous-shock-layer  (VSL3DNQ)  code  has  been 
modified  for  analyzing  the  low-density,  nonequilibrium  flow  past  hypervelocity 
vehicles.  This  code  is  robust,  dependable  and  fairly  accurate. 

2.  Recently  obtained  surface  and  shock-slip  boundary  conditions  are 
implemented  to  account  for  the  low-density  effects. 

3.  The  governing  equations  are  solved  in  the  computational  domain  with  a 
uniform  grid.  The  grid-transformation  employed  gives  finer  resolution  at  both 
surface  and  shock  and  does  not  add  any  numerical  dissipation  to  the  results. 

4.  Input  shock  shape  obtained  from  the  viscous  nonequilibrium  flow 
results  in  a  superior  convergence  history,  especially  for  the  low-density 
flows. 

5.  Very  good  agreement  is  obtained  with  the  deary's  perfect-gas  data 
for  the  high  Reynolds  number  case. 

6.  Results  for  the  zero-deg  angle-of-attack  case  are  in  excellent 

agreement  with  those  obtained  from  an  earlier  developed  2D  nonequilibrium 

code . 

7.  Excellent  agreement  is  obtained  with  the  RAM  C  flight  data  for  the 
low-density  high-energy  flow  conditions. 
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3. 


Present  results  compare  favorably  with  those  obtained  from  the 


Direct-Simulation  Monte  Carlo  (DSMC)  method  for  the  low-density  hypersonic 
f lows . 

Finally,  despite  its  certain  limitations,  viscous  shock-layer  method  has 
been  shown  to  be  computationally  efficient  and  a  valuable  tool  for  analyzing 
the  low-density  flows.  With  the  method  suggested  here  for  obtaining  the  input 
shock  shape  for  the  3D  nonequilibrium  viscous  flows,  the  VSL  method  becomes 
self-starting  unlike  the  parabolized  Navier-Stokes  (PNS)  methods.  The  VSL 
method  has  been  used  for  obtaining  the  starting  profiles  for  some  of  the  PNS 
calculations.  The  surface  slip  conditions  presented  here  may  also  be  used  for 
low-density  calculations  employing  the  Navier-Stokes  equations.  Further,  if  a 
shock-fitting  approach  is  used  for  such  calculations  to  keep  the  computer 
storage  requirement  reasonably  small,  the  present  shock-slip  formulation  can 
be  utilized. 
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Figure  1  BODY  AND  SHOCK- 
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I'igure  2a.  SHOCK- WAVE  NORMAL  AND  TANGENT  VECTORS  AND  ROTATION  ANGLES  ABOUT  N  AND 
e*  VECTORS 
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Figure  2c.  VELOCITY  COMPONENTS  AND  ANGLES  EMPLOYED  FOR  SHOCK-BOUNDARY  CONDITONS 
IN  A  MERIDIONAL  =  CONSTANT)  PLANE 
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Figure  3.  SOLUTION  OF  VISCOUS  SHOCK  LAYER  EQUATIONS  ON  THE  BODY  (  OR  CONICAL  PART  ) 
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Figure  4.  SOLUTION  OF  NOSE  REGION  (  OR  SPHERICAL  PART  ) 
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Figure  5.  COORDINATES  AND  BODIES  AT  ZERO  AND  a-DEG  ANGLES  OF  ATTACK 
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COMPARISON  OF  SURFACE  HEAT  TRANSFER  RATIO  BETWEEN  EXPERIMENTAL 
AND  NUMERICAL  RESULTS  FOR  A  15“  SPHERE-CONE 
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Figure  7.  COMPARISON  OF  SURFACE  HEAT  TRANSFER  AT  DIFFERNT  CROSSFLOW  PLANES 
FOR  A  15°  SPHERE-CONE  WITH  AND  WITHOUT  SLIP 


DIFFERENT  APPROACHES 


MONEQUILIBRIUM  FLOW 


OUTPUT  SHOCK  LAYER  THICKNESS  FOR  A 


Figure  10  .  SMOCK  CONVERGENCE  HISTORY  FOR  A  9°  SPHERE-CONE  USING 
VISCOUS  NONEQUILIBRIUM  INPUT  SHOCK  SHAPE 
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INVISCID  EQUILIBRIUM  INPUT  SHOCK  SHAPE 


Figure  12  .  COMPARISON  OF  HEAT  TRANSFER  COEFFICIENT  FOR  A  9°  SPHERE- 
CONE  WITH  TWO  DIFFERENT  INPUT  SHOCK  SHAPES 
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COMPARISON  OF  VELOCITY  DISTRIBUTION  AT  s=0.0  FOR  A  20° 
SPHERE -CONE  BETWEEN  2-D  AND  3-D  NONEQUILIBRIUM  CODES 
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Figure  14  .  COMPARISON  OF  TEMPERATURE  DISTRIBUTION  AT  6=0.0  FOR  A  20' 
SPHERE-CONE  BETWEEN  2-D  AND  3-D  NONEQUILIBRIUM  CODES 
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Figure  15  .  COMPARISON  OF  MASS  FRACTION  DISTRIBUTION  OF  ATOMIC  OXYGEN  AT 
s=00  FOR  A  20*  SPHERE-CONE  BETWEEN  2-D  AND  3-D  MOH- 
EQllll  IBRIUM  CODES 
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Figure  16  .  COMPARISON  OF  MASS  FRACTION  DISTRIBUTION  OF  MOLECULAR  OXYGEN  AT 
s=O.Q  FOR  A  20°  SPHERE-CONE  BETWEEN  2-0  AND  3-0  NON- 
EQUILIBRIUM  CODES 


Figure  17  .  COMPARISON  OF  MASS  FRACTION  DISTRIBUTION  OF  ATOMIC  NITROGEN  AT 
s=Q.O  FOR  A  20°  SPHERE-CONE  BETWEEN  2-0  AND  3-D  NON- 
EGUIUBRIIJM  CODES 
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.  COMPARISON  OF  MASS  FRACTION  DISTRIBUTION  OF  MOLECULAR  NITROGEN  AT 
s=0.0  FOR  A  20®  SPHERE -CONE  BETWEEN  2-D  AND  3-D  NON- 
EQUIIIBRIIJM  COOES 
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Figure  20  .  SURFACE  HEAL  TRANSFER  COEFFICIENT  FOR  A  20*  SPIIERE-CONE  AT 
DIFFERENT  CROSSFLOW  PLANES  -  MO  SI  IP 


MASS  FRACTION  DISTRIBUTION  FOR  ATOMIC  OXYGEN  AT  s=2.02  FOR  A 
20°  SPHERE-CONE  AT  DIFFERENT  CROSSFLOW  PLANES-NO  SLIP 
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MASS  FRACTION  DISTRIBUTION  FOR  ATOMIC  NITROGEN  AT  s=2.02  FOR 
20°  SPHERE-CONE  AT  DIFFERENT  CROSSFLOW  PLANES-NO  SLIP 
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FRACTION  DISTRIBUTION  FOR  MOLECULAR  NITROGEN  AT  S=2.02  FOR 
1  SPHERE -CONE  AT  DIFFERENT  CROSSFLOW  PLANES-NO  SUP 


ELECTRON  HUMBER  DENSITIES  AT  s=2.02  FOR  A  20'  SPHERE-CONE 
AT  DIFFERENT  CROSSFLOW  PLANES 


COMPARISON  OF  UFA!  TRANSFER  COEFFICIFNT  FOR  A  9 
SPHERE-CONE  AT  DIFFERENT  ANGLES  OF  ATTACK 
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Figure  JO.  COMPARISON  OK  SURFACE  HEAT  TRANSFER  COEFFICIENT  AT  DIFFERENT 
ALTITUDES  FOR  A  5°  BLUNTED  CONE  WITH  DIRECT  SIMULATION 
MONTE  CARLO  MZTHOD 
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altitude:,  Km 

Figure  31.  COMPARISON  OF  SURFACE  PRESSURE  COEFFICIENT  AT  DIFFERENT 
ALTITUDES  FOR  A  5°  BLUNTED  CONE  WITH  DIRECT  SIMULATION 
MONTE  CARLO  METHOD 
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32.  COMPARISON  OF  STAGNATION  REGION  TEMPERATURE  PROFILES  (WITH  SLIP) 
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Figure  33.  COMPARISON  OF  SUFFACE  HEAT  TRANSFER  DISTRIBUTION  FOR  A  9' 
SPHERE-CONE  (WITH  SUP) 
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Figure  3^.  COMPARISON  OF  SURFACE  HEAT  TRANSFER  COEFFICIENT  FOR  A  9' 
SPHERE-CONE  (WITH  SIJP) 


Tangency  Points 


COMPARISON  OF  SHOCK  LAYER  THICKNESS  GRADIENT  FOR  A  9' 


figure  36  .  COMPARISON  OF  SMM-FRlCT)ON  COEFFICIENT  FOR  A  9 
SPHERE -CONE  WITH  AND  WITHOUT  SLIP 


COMPARISON  OF  HEAT  TRANSFER  COEFFICIENT  FOR  A  9 


C£ 

a 

a: 

a 

<j) 

n 

CO 

.C 

o 

d 

n 

<N 

o 

a 

r- 

K 

in 

o 

<0 

GO 

ai 

CN 

m 

m 

i— 

o 

II 

II 

a 

il 

ii 

ii 

d 

O  t  .  ■  .  *  • 9  || 

2  2  l—  3  Q£  h—  u 


Figure  38  .  COMPARISON  OF  SURFACE  PRESSURE  COEFFICIENT  FOR  A  9' 
SPHERE-CONE  WITH  AND  WITHOUT  SUP 
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Figure  39.  COMPARISON  OF  STREAM  WISE  VELOCITY  DISTRIBUTION  FOR  A  9 
SPHERE-CONE  AT  a=8.8  WITH  AND  WITHOUT  SLIP  FOR  <t>=0° 
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Figure  40.  COMPARISON  OF  TEMPERATURE  DISTRIBUTION  FOR  A  9°  SPHERE-CONE 
AT  6=8.8  WITH  AND  WITHOUT  SUP  FOR  <t>=0° 
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COMPARISON  OF  DENSITY  DISTRIBUTION  FOR  A  9°  SPHERE-CONE 


Figure  42.  COMPARISON  OF  MASS  FRACTION  DISTRIBUTION  FOR  ATOMIC  OXYGEN  FOR 
A  9°  SPHERE-CONE  AT  e=8.8  WITH  AND  WITHOUT  SLIP  FOR  4^0 
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COMPARISON  OF  MASS  FRACTION  DISTRIBUTION  FOR  MOLECULAR  OXYGEN  FOR 
A  9°  SPHERE-CONE  AT  a=8.8  WITH  AND  WITHOUT  SUP  FOR  4»=0 
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COMPARISON  OF  MASS  FRACTION  DISTRIBUTION  FOR  ATOMIC  NITROGEN  FOR 
A  9“  SPHERE -CONE  AT  a=8.8  WITH  AND  WITHOUT  SUP  FOR  4-=0 


0.79 


-86- 


COMPARISON  OF  MASS  FRACTION  DISTRIBUTION  FOR  MOLECULAR  NITROGEN 
A  9*  SPHERE-CONE  AT  6=8.8  WITH  AND  WITHOUT  SLIP  FOR  <T=0 
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Figufu  T6.  COMPARISON  OF  MASS  FRACTION  DISTRIBUTION  FOR  NITROGEN  OXIDE  FOR  A 
9°  SPHERE-CONE  AT  e-8.8  WITH  AND  WITHOUT  SUP  FOR  T>=0 


